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Abstract. A partial description of Earth’s broad scale, core-source magnetic field has been developed 
and tested three ways. The description features an expected, or mean, spatial magnetic power spectrum that 
is approximately inversely proportional to horizontal wavenumber atop Earth’s core (Stevenson, 1983; 
McLeod. 1985; 1996). This multipole spectrum describes a magnetic energy range; it is not steep enough 
for Gubbins (1975) magnetic dissipation range. Temporal variations of core multipole powers about mean 
values are to be expected and are described statistically, via trial probability distribution functions, instead 
of deterministically, via trial solution of closed transport equations. The distributions considered here are 
closed and neither require nor prohibit magnetic isotropy. The description is therefore applicable to, and 
tested against, both dipole and low degree non-dipole fields. 

In Part I, a physical basis for an expectation spectrum is developed and checked. The description is 
then combined with main field models of twentieth century satellite and surface geomagnetic field 
measurements to make testable predictions of the radius of Earth’s core. TTie predicted core radius is 0.7% 
above the 3480 km seismologic value. Partial descriptions of other planetary dipole fields are noted. 

1 INTRODUCTION & BACKGROUND 

The primary physical source of the geomagnetic field is widely held to be electric current flowing in 
Earth s electrically conducting, ferro-metallic liquid outer core and solid inner core. Other sources include 
weaker current in the resistive, ferro-magnesian silicate and oxide mantle, magnetization in the colder crust, 
electric current in the saline hydrosphere, and external currents in the ionosphere above and magnetosphere 
beyond. Solenoidal magnetic induction B within the source regions can be represented mathematically via 
poloidal and toroidal scalars (see, e.g.. Backus, 1986). On and above Earth’s roughly spheroidal surface, 
however, the field of deep internal origin can be represented by the negative gradient of a scalar potential V 
that satifies Laplace’s equation (B = -VV, V^V = 0). At time t and position r in geocentric spherical polar 
coordinates (radius r, colatitude 0, and east longitude <|)), the spherical harmonic expansion of zero mean V 
is well-known to be 

" n+1 " 

V(r,t) = a E (a/r) I [gn"^(t)cosm<t) + h^f^COsinmitilP f"(cos0) , n) 

n=l m=0 ^ ^ 

where Pj^ is the Schmidt-normalized associated Legendre function of degree n and order m, and [gf|^(t), 

h„"’(t)] are Schmidt-normalized Gauss coefficients at reference radius a, here taken to be 6371 .2 km. 

Gauss coefficients have long been estimated by a weighted least squares fit to measured geomagnetic 
data (see, e.g., Langel, 1987). This usually requires truncating the sum over degrees at finite N. Accurate 
statistical information about Gauss coefficients could supplement such regularizing truncation, increase the 
accuracy of estimated coefficients, and improve the reliability of associated uncertainty estimates (Gubbins, 
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1983; McLeod, 1986; Backus, 1988). Lacking sufficient prior information, a statistical hypothesis may be 
advanced and tested against observations. The physical significance of such tests depends upon the 
physical basis of the hypothesis tested. Discrepancies between theory and observation often lead to more 
accurate hypotheses and measurements. So let us consider such hypotheses, develop their theoretical 
predictions, and test the predictions against geomagnetic and paleomagnetic observations. 

Many hypotheses concern the internal spatial geomagnetic power spectrum, denoted via Rj^. Lowes 
(1966) and others (see Cain et al., 1989) show the mean square potential field represented by spherical 
harmonics of degree n, averaged over a spherical shell of radius r enclosing the sources, to be 

2n+4 ^2 2 

R„(r,t) = (n+l)(a/r) I ). (2) 

m=0 

For vacuum magnetic permeability the magnetic energy per unit shell thickness due to an equivalent, 
geocentric, n^ degree multipole source equals 27tr^Rj^/pQ, The latter form a discrete spectrum with units 
of spatial power or force. "Multipole powers" R^^ have SI units of (Tesla)^ and collectively sum to the 
mean square field on the shell. Invariant under coordinate rotations, but not translations, the Rj^ offer 
isotropic, inhomogenous measures of an evidently anisotropic and heterogenous field. The spatial power 
spectrum of the secular variation, denoted Fj^(r,t), is obtained by replacing the Gauss coefficients in (2) with 
secular variation (SV) coefficients (3jgj^^,5^hj^^); note d^^R^ ^ ^ ^t^n* 

Multipole powers R^ (a, 1980) computed from coefficients fitted to Magsat data are interpreted in 
terms of a predominantly core-source spectrum R^^^ for n ^ 12 and a predominantly crustal-source 
spectrum R^^^^ for n > 16 (Langel & Estes, 1982). Corrections for lithospheric magnetization and finite 
mantle resistivity are arguably small only at low degrees and long periods, respectively; therefore, to ease 
description of a field originating in a core of mean radius c, use of (2) is here limited to n < 12, super- 
annual SV, and r>c. Seismologic estimates of c vary, but are denoted Cg = 3480km (Dziewonski & 
Anderson, 1981; Kennet, Engdahl & Buland, 1995). 

Several functional forms have been suggested for R^^^. The theoretical core-source geomagnetic 
spectrum of McLeod (1985; 1996 equation 20a) earns special attention, not only becasue it can be closely 
fitted to observational but because it was obtained from consideration of advective and diffusive effects 
on a core field. McLeod’s rule states that the spatial magnetic power spectrum of the core-source field, 
Rnc(a)» is approximately proportional to (n + 1/2)' ^(c/a)^" for finite degrees 2 < n < N’. Here we advance 
and test the hypothesis that this form approximates geologically long time average (over ~50 Myr), or 
expectation, multipole powers throughout a finite magnetic energy range 1 < n < Ng. Evidently, 
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Ng ^ 12. Natural fluctuations about such mean values are to be expected from time-dependent magnetic, 
mass, momentum and energy transport within Earth’s core. Over longer time intervals, changes in outer 
core boundary conditions, composition, and stratification may change itself. 

In Part I, after summarizing some previous work with the spectrum, a simple physical model of the 
core-mantle boundary region is used to develop an expected core multipole spectrum. The model features 
laterally decorrelated magnetic transport beneath a thin viscous sub-layer. This specifies the form of SV 
spectrum A relation between effective radial wavenumbers and horizontal wavenumbers in the 

vertical field component is needed to specify {RncJ* One such relation is advanced and tested against 
observational spectral ratios with satisfactory results. It yields an expected spectrum similar to McLeod’s. 
A trial distribution function for is advanced and discussed, mainly to guide the paleomagnetic 
predictions and tests presented in Part II. Geomagnetic predictions of the mean radius of Earth’s core based 
on McLeod’s rule are then tested against seismologic values and found to be quite accurate. 

2 SOME PREVIOUS SPECTRAL FORMS 

Exponential Spectra. The geometric attenuation factor, in (2), is a power law in radius; 

yet it is exponential in degree. Following Lowes (1974), Rjj(a) has often been approximated by the 
exponential A*(c*/a)^^'^. Linear regression through observational values of Iog[Rj^(a)] gives 

Iog[R*(a)] = [21og(c*/a)]n + [log(A*) + 41og(c*/a)] . (3) 

The slope of this line implies the radius c* at which downwardly continued spectrum R*(a)[a/r]^*^'*'^, 
denoted Rp(r), would appear independent of n. A graph of log[Rf^(r)] as a function of n is a straight line 
that "levels off at r = c . If Rfj(r) were extrapolated to arbitrarily high degree, then its sum over n would 
diverge at r < c*; therefore, such extrapolation is physically invalid at r < c*. This might be due to failure 
of the potential field representation, so c* might be the minimum radius of a sphere containing the sources. 
Lowes’ (1974) spectrum, fitted to degrees 1 through 8, levels off about 480 km below c^. Modem dipole 
power Rj exceeds Rj as well as higher multipole powers; it has often been excluded from subsequent 
regressions, which raises c* towards c^. 

Langel & Estes (1982) interpret Rj^(a) from their degree 23 model MGST 10/81 of Magsat data in 

terms of dominant core dipole, core non-dipole (2 < n < 12), and crustal fields (n > 16). Their non-dipole 

♦ ■ ♦ 
core spectrum levels off 174 km below c^. Their crustal spectrum Rj^^ levels off 83 km below a. 

Voorhies (1984) used only the mean square radial field component per harmonic degree from model 

MGST 10/81, which is (n-»-l)Rj^/(2n-i-I). His non-dipole core spectmm B*^^(n) levels off 200 km below Cg 
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and was used to argue against narrow scale* intensely magnetized core spots at the core-mantle boundary. 

His crustal spectrum B*^^(n) levels off 9.9 km below a* perhaps suggesting crustal rather than deep 

lithospheric sources. About 10% of Rj 2 was considered crustal in origin. 

Cain et al. (1989) used their degree 63 numerical integration model M07AV6 of Magsat data to 

obtain a non-dipole core spectrum that levels off 73 km below c^; a crustal spectrum that flattens out 14 km 

below the reference sphere; and a satellite altitude noise level. 

Constable & Parker (1988) chose a non-dipole spectrum which is exactly flat at the core surface to 

guide construction of a statistical model of Gauss coefficients. Except for the manifestly anisotropic dipole* 

and perhaps quadrupole, normal probability distribution functions (PDFs) with m-independent variances 

were assigned to the coefficients. The n-dependent variances were obtained from a level spectrum at Cg 

fitted to observational of model GSFC 9/80 (Langel* et al., 1982). The model passed a Kolmogorov- 

Smimoff test for degrees 2-8, 
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Hulot, LeMouel & Wahr (1992) inferred a non-dipole core spectrum that levels off about 235 km 

If 

below Cg. The stationary isotropic statistical model of Hulot & LeMouel (1994) also assigns normal PDFs 
with m-independent variances to non-dipole Gauss coefficients. Variances obtained from decrease 
quite rapidly with n due to the depth of c* below Cg. They also offer impressive non-dipole Kolmogorov- 
Smimoff test results, note that isotropic models give chi-squared PDFs for normalized non-dipole powers, 
and further suggest a PDF for the axial dipole coefficient that combines two Gaussian distributions: one for 
normal and one for reversed polarities (see McFadden & McElhinny, 1982). 

The statistical models cited above allow for the dominant geometric attenuation of R^qW with radius 
and seem to offer unobjectionable, albeit different, statistical hypotheses. Yet physical reasons to expect 
non-interdependent Gauss coefficients, isotropy, equal degree variances at some radius c*, and a purely 
exponential spectrum - all while granting special exemptions to the axial dipole - are elusive. 

For example, with magnetic permeability \x and electric conductivity o simply treated as if uniform in 
a spherical fluid outer core moving at velocity v(r,t), analysis of the well-known induction equation, 

3(B= Vx(vxB) + (nor^V^B , (4) 

indeed shows that magnetic diffusion of different spherical harmonic modes proceeds independently; yet 
poloidal and toroidal magnetic modes are cross-coupled by heterogeneous fluid motion - the motion 
needed to maintain a core-source field against resistive decay of its source currents over geologic time via 
dynamo action (see, e.g., Larmor, 1919; Elsasser, 1946; Gubbins & Roberts, 1987). Moreover, Earth’s 
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present dipole is evidently anisotropic and has power far greater than obtained by extrapolating the 
foregoing exponential spectra to the first degree. Furthermore, no firmly established, fully developed 
statistical theory of conceivably turbulent core magnetohydrodynamic transport appears to demand equal 
mean magnetic pressure per harmonic degree at radius c (level )) instead of, say, equal total 
magnetic energy outside the core per harmonic degree and order (level (n+l)(2n+l)Rp(c)), equal mean 
radial gradient in magnetic pressure per harmonic degree atop the core (level R^(c)/(n+2)), or some other 
spectral form at degrees below 13. Alternatives to a purely exponential spectrum, and to statistical 
presumption of isotropy, independence, and normally distributed coefficients, are available and have been 
used to describe geo-paleomagnetic field behavior. 

Power Law Spectra by the Core-Mantle Boundary. The exponential decay of Rj^(r > > c ) with 

n is faster than needed for finite mean magnetic energy density (faster than n'^"^ with e > 0). Indeed, the 
finitude of Ohmic dissipation within the core and Gubbins’ (1975) famous expression for the minimum 
value thereof imply that, for degrees in a magnetic dissipation range defined by n > Nj^, 

Rn^(a) < Kq n-2-8 (c/a)^" , (5) 

where Kq is a constant and 5 > 0. Though finite, Nj^ might be very much greater than 12, 

Much work on idealized hydromagnetic turbulence has considered the magnetic energy spectrum 
E(k) as a function of Cartesian wavenumber k (see, e.g., Gubbins, 1974; Pouquet, Frisch & L^orat, 1976; 
Moffat, 1978; Krause & Radler, 1980). Although the proportionality E(k) a k'^^^ is indicated for the 
inertial subrange of three-dimensional, homogeneous, isotropic, incompressible hydromagnetic turbulence, 
Pouquet et al. (1976) found that kinetic helicity injection leads to an inverse cascade of magnetic helicity 
and E(k) a k*^ at low k. Citing these findings, Stevenson (1983) suggested that a power law spectrum is 
possible - albeit difficult to apply due to anisotropy. His arguments indicate k‘^ and he noted some 
similiarity between an n"^ spectrum and observational R^^ downwardly continued to radius 0.55a for n < 8. 

Define horizontal wavenumber k^j^ s [n(n-i-l)/c^] via the surface Laplacian. If E(k)ak'^, if 

E(k) a {Rnc(c)}’ and if k^ a j^en {Rnc(c)l « « [n(n+l)]-^^^. The latter form is closer to 

McLeod’s (n + 1/2)’^ than Stevenson’s n"^ relation, notably at degrees 1 and 2. 

Gubbins & Bloxham (1985) discuss solution norms for supplementary regularization of core-source 
field models. Each such norm corresponds to a downwardly continued spectrum that is inversely 
proportional to a polynomial of n. Norms compatible with (5) are needed for n > Nq. Solution norms may 
also be used to construct a smooth, or damped, core field models for n < Nj) (see also Shure, Parker & 



Backus (1982); Langel (1987); and Backus (1988)). Such models have many uses; however, a desire fora 
smooth model ought not be mistaken for measured data. To avoid built-in assumptions about c or 
smoothed field models are not used in this paper. 

McLeod’s (1985; 1994, 1996) theoretical core-source spectrum for degrees 2 < n < N’ is written 
^nc^®^ (c/a)^" , 

where = 5x 10 ^ nT2 The downwardly continued spectnim R^(t) is proportional to (n + 1/2)- 1 at the 
core surface; there is no radius at which it levels off. McLeod (1996) derived ( 6 ) from functional forms of 
(i) the core-source S V spectrum and (ii) the temporal geomagnetic power spectrum P„(co) that depends 
upon temporal frequency to as well as harmonic degree. The form of F„, is akin to that caused by random 
lateral advection of magnetic field lines at the top of a high conductivity liquid core. The form of P^fco) is 
appropriate to a two time-scale model of the two processes, magnetic flux diffusion and fluid motiol. that 
change the core field according to the induction equation ( 4 ). 

Voorhies & Conrad (1996) questioned the existence of physical sources giving a level spectrum, the 
exclusion of dipole power from estimates of rJ,, and the idea that < 14. ExaminaUon of the spectrum 
from the degree 60 Magsat model M102189 of Cain. Holter & Sandee (1990) confinned that dipole power 
is large, and quadrupole power is small, compared with a trend at degrees 3 - 12 . TTiis was thought to be an 
artifact of geologic undersampling: values for Rj and R 2 at 1980. being calculated from but 3 and 5 Gauss 
coefficients respectively, could be fairly far from multi-million year mean values. If so, perhaps a key to 
past dipole behavior can be found in present higher degree multipoles. We tried the power law form at c^ 

= Ky 3 ‘‘"ear function of ln(n) to InfR^fc^)] for degrees 3-12. We found 7 to be 

0.94. Crustal sources may lead to a slight underestimate of y, so we set y to 1, estimated Ky alone, and 
found the summed squared residuals per degree of freedom to be less than when y and Ky were coestimated. 
Moreover, the extrapolation to dipole power appeared much more satisfactory than for a plain exponential 
form. With neither theoretical nor experimental reasons to assume a level spectrum at either c^ or c*. and 
seeing fair predictions for R j and R 2 , we learned to expect a core-source spectrum of the form 

RncCrtO = {R„l(r)} = Kj (c/r)2n-H4 

for 1 < n < Ng, where the curly brackets represent expectation value and Ng is the finite, maximum degree 
of the magnetic energy range (1 < Ne < N^). Whether or not a Stevenson (1983) relation like ( 7 ) describes 
t’ Myr means, it does not specify changes in c, y, or Ky due to planetary evolution over several Gyr. 


If (7) accurately describes the time-averaged core-source field, then natural deviations [R„(;(t) - 
(Rnc)] eventually relax to zero and perhaps change sign. Finding Rj somewhat high, and R 2 low, at 1980 
led us to predict relaxation of these deviations. We tested this against the main field and SV coefficients of 
undamped model GSFC 12/83 (Langel & Estes, 1985). Rj is decreasing and R 2 is increasing as predicted. 
Moreover, for all orders m, Ojgi"’)/gi"’ < 0 and Otg 2 '")/g 2 "’ > 0 for both g^"* and h„"™. The plain 
chance of such perfect correlation of signs is 1/8 for the dipole, 1/32 for the quadrupole, and 1/256 for all 
eight coefflcients together. This otherwise remarkable coincidence can be viewed as efficient relaxation of 
the core-source field towards geologic mean values like (6) or (7). 

3 STATISTICAL PHYSICS OF THE CORE-SOURCE FIELD 

For simplicity, the transition from core to mantle is approximated by a sharp material interface: a 
jump discontinuity in molecular material properties with radius. Both hydrodynamic and macroscopic 
magnetic stress tensors are, in contrast, treated as continuous across the interface to ensure finite divergence 
of the magnetohydrodynamic stress and finite viscous and Lorentz force densities. If the fluid wets the 
mantle and is of approximately Newtonian rheology, then motion of the viscous fluid vanishes at the 
interface. The molecular kinematic shear viscosity of the liquid metal v is very small, about 3x 10"^ m^/s 
(Poirer, 1988; also see Lumb & Alldridge, 1991), so lateral motion ought not be neglected beneath a thin 
viscous sub-layer. A laminar sub-layer as thin as the Ekman depth dg or the Hartmann depth djj of about 8 
cm or 20 cm, respectively, need not imply that diffusion is unimportant in the sub-adjacent region, which 
might be a thicker boundary layer that has not usually been distinguished from an underlying main stream. 

Traditional analyses yield a thin, weak visco-magnetic boundary layer that neither absorbs 
appreciable normal fluid flow nor generates much electric current (Roberts & Scott, 1965; Backus, 1968; 
Ball, Kahle & Vestine 1968; Hide & Stewartson, 1972; Benton, 1981; Gubbins & Roberts, 1987; Benton, 
1992). Advection of the magnetic field at the top of the main stream is thus mainly lateral; moreover, the 
jump in magnetic field across such a layer is so small as to be negligible for the radial component Bj.. The 
radial component of (4) has thus been used to interpret SV in terms of broad scale fluid flow, and 
occasionally flux diffusion, at the top of the main stream. Constrained inversions of modem geomagnetic 
secular change indicate a main stream speed U of about 7.5 km/yr (see, e.g., Voorhies, 1995). Such 
inversions often damp out poorly resolved narrow scale motions. Yet the corresponding main Reynolds 
number Uc^/v of about 3x 10^ is so much greater than the boundary number Udg/v of about 60 as to 
suggest some ephemeral small scale motions occur in the core. Secular change caused by such motions 
might be described via a pseudo random walk of magnetic field line footpoints. 




For example, Benton’s (1992) lateral flow scale of 4.8 km might describe the width of fronts between 
broad regions of more uniform core flow, but it might also describe a seething mass of short-lived, 
rotationally polarized hydromagnetic eddies. Such eddies would individually induce unobservably narrow 
scale field variations (degrees of roughly 3,000), yet may collectively contribute to observable broad scale 
SV. Statistical parameterization of this contribution could be an efficient, useful, and informative 
complement to deterministic modeling. 

There is a direct statistical approach to lateral magnetic transport at the base of the viscous sub-layer: 
the equivalent source approach of McLeod (1996). This yields a kinematically unbiased ensemble mean 
SV spectrum (Fj^) (section 3.1). Whether or not main stream dynamics cause geologic time averages to 
differ from such an ensemble average, observable core signal diffuses through the sub-layer and resistive 
mantle. To do so, magnetic change originating deeper in the core, be it induced by eddies 2 or 22,000 km 
across, induces electric current, hence excess field line curvature (V^B/I B I ), in the sub-layer (section 3.2). 
One hypothesis about such curvature is advanced and tested by observation and analysis of the emergent 
core field (section 3.3). It yields an expected spatial power spectrum for the broad-scale core-source field 
{Rncl similar to McLeod’s rule. It does not specify spectral variance {[R^^ - higher 

moments. Closure is sought via trial PDFs for Rj^^ that (i) have variances determined by the mean and (ii) 
neither require nor prohibit isotropy (section 3.4). 

3.1 Statistical Kinematics of the Secular Variation. 

Following McLeod (1996), quasi-static lateral magnetic transport causes differential exterior secular 
change AB equivalent to differential dipole moments Ad^ scattered atop the main stream (r = c"). To see 
this, recall that a single magnetic flux vector B^dA at fixed position on the surface c~ of the source 
region acts as the point source of a dipole field with moment proportional to B^dA. The magnetostatic 
field at position x due to this equivalent source is well-known (see, e.g., Jackson, 1975 equation 5.64). 
Infinitesimal quasi-static lateral displacement Ax of this single magnetic vector, with no change in 
orientation and magnitude, causes a net change in the exterior field equivalent to a differential quadrupole 
moment at x^ + Ax/2. More generally, there is a magnetic vector at each position on the source-surface; 
lateral transport can replace the vector at x^ with an adjacent vector of slightly different orientation and 
magnitude; and the change in the exterior field is equivalent to that of a differential dipole moment Ad^ at 
x^. Given many equivalent source changes Adj at x- on c’, the total change in the exterior field AB at 1x1 
> c' follows by superposition. In the continuum, differential surface magnetic moment density at x’ on c' 
replaces discrete Adj as the equivalent source of exterior secular change AB(x,t), which follows by 
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integration over x\ Random lateral transport of the core-source field thus produces a corresponding change 
in the field above. 

Elements of the dyad formed by two differential dipole moments at two well separated points 
[Adj][Ad 2 ]^ may be either positive or negative. The average over a kinematically unbiased ensemble of 
dyads gives zero cross-correlations, but non-zero auto-correlations. Then the expected spatial power 
spectrum for broad-scale secular variation is equivalent to the spectrum for laterally uncorrelated, 
randomly varying dipole moments on the shell of radius c’ 

{F„j.(r > C-) ) = C n (n + 1/2) (n + 1) (c-/r)2n+^ . (8) 

This differs slightly from McLeod (1996, equation (11)) because changes in horizontal as well as radial 
components of the equivalent core-source field can contribute to Fj^^. (Consider rotation about the vertical 
of a horizontal equivalent source on the equator). Laterally decorrelated point sources may offer the 
roughest physical model of SV, yet the sum of geometrically attenuated cubic spectrum (8) over n 
converges on all spheres above the equivalent source layer. 

When (8) is used to describe a geologic time average, natural fluctuations - {Fp}] are described 
as transient lateral cross-correlations between equivalent SV sources. Failure of (8) at spatial scales 
corresponding to n > N£ may in part be due to persistent mesoscale cross-correlations. Moreover, 
description of non-stationarity via (8) is limited to changes in amplitude C and in c’ between geologic time 
averaging intervals. Because C is proportional to mean square differential dipole moment {(Ad| /At)^}, it 
should, by (4), tend to increase with mean transport speed and mean field intensity. 

Remarks. The quasi-magnetostatic derivation of (8) is so fundamental that it does not require use of 
Faraday’s induction equation. Magnetic transport equation (4) obtains when (i) Galilean invariance and 
Ohm’s law, hence transport of non-uniform field B(r,t) by non-uniform fluid velocity v(r,t), replace 
transport of equivalent sources; (ii) Ampere currents of density J are the true core-sources; (iii) Maxwell 
displacement currents are filtered out via the quasi-steady approximation (V«J = 0); and (iv) p and a are 
treated as if uniform scalars. For stationary reference mass density PQ(r), the anelastic approximation of 
Gubbins & Roberts (1987) is replaced with V^p^v = 0. Then expansion of a rising fluid parcel (V«v = 
-v*VlnpQ > 0) tends to reduce its magnetic flux density and non-stationarity is discussed separately. At the 
base of the sub-layer, omitting terms proportional to small Vj. (hence V«v), the radial component of (4) is 


9jBr+ v.VjB, = 


Bj.a^v^ + ()iac)'^v2rBj., 


( 9 ) 


where denotes the surface gradient oj>erator. 

Observable changes in the core-source field caused by advection and by downwelling 

(BrdfVr) are indistinguishable from quasi-static lateral transport and confluence of equivalent sources. The 
broad scale part of these changes are here described via (8). Ohmic resistance to Ampere currents appears 
as magnetic diffusion in (9). Observable changes in the field due to diffusion are indistinguishable from 
changes in equivalent source moments. Whether or not magnetic diffusion alters (8), as seems likely for n 
^ and during laterally correlated free decay, it is essential in the sub-layer where v falls to zero. 

The lateral decorrelation used to obtain (8) by no means prohibits turbulent dynamo action at depth; 
indeed, the more precisely correlated fluctuations balance diffusion, the better (8) should describe any 
residual SV. Although changes in horizontal field components also contribute to (8), the normal field 
component on the two boundaries (r = c, r ^ still detennines the exterior potential field. So long as Bj. 
remains effectively continuous across the sub-layer, so does then jumps in Vg®Bj.v and (jiGr)'^ V^Bj. 
across the sub-layer are the same by (9). The statistical description of broad-scale SV (8) is, by continuity, 
therefore applied atop the core-mantle interface. 

3.2 Apparent Scale Heights and Characteristic Time-Constants 

At the top of the viscous sub-layer, v vanishes by the no-slip boundary condition. There (4) reduces 
to the magnetic diffusion equation with radial component 

a^r = (4<Jc)-lv2rBr = . (10a) 

'y ^ 

Here Bj. denotes the jump in 3j. Bj. across the core-mantle interface. Clearly Bj.* is the core-source 

correction (Bj. + a^V) for r < c. It vanishes at c by continuity. From (10a) alone, B^* would seem 

comparable in magnitude to Bj.(c) at a depth below c, denoted kj.'^ of about [<[Bj.(c)]2>/<[paaj 
2 1/2 

BfCc)] >] ' , where the angle brackets denote the average over the sphere. This is an apparent scale height 
because (10a) does not hold in moving fluid below c. At Magsat epoch 1980, root mean square (rms) 
values for the downwardly continued broad scale radial main field and SV are about 3 gauss and 3 pT/yr; 
with n = Hq, assuming o = 5x10^ S/m would return an apparent scale height of about 70 km. 

Denote real, Schmidt-normalized spherical harmonic coefficients of Bj.* within the core by [Gp'^(r,t), 
Hn^’fr.t)]. Omitting corrections for mantle conductivity and crustal sources, (10a) implies 

(n -h l)(a/c)"+2a^g^m ^ (^^)-l 3 ^ 2 G^m 
(n -t- l)(a/c)"+2a^h„f" = (pa)-l a^2H^m 
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(This plausibility argument neither suggests nor requires a weak field deep in the core). More generally, 
when the physical dependencies of [kj.|^(t)]^ on underlying hydromagnetic interactions, and on magnetic 
modes contributing to many other multipole powers Rjj,(c,t) for j * n, dominate its dependence on Rpj,(c,t), 
then (13) can be a useful approximation. 

Statistically, for any set of PDFs assigned to (£„*", hn™) and for any set of PDFs assigned to 
f^r^^n"*’ ^r^^n'^J at r = c, define the ratio of expectation values {Bj.p^(n)}/{[9|.^ Bj.*(n)]^} to be kj.^"^; 
kj^ is the ^ective radial wavenumber for the radial field configured in harmonics of degree n. As in (1 la- 
d), one obtains 

(Fnc) = (pok„2)2 

To the extent that expectation values correspond to geologic time averages* (14a) corresponds to (12) rather 
than (13); indeed* only during intervals when correlated fluctuations are negligible compared with the 
product of means could carefully averaged apparent scale heights {[kj^(t)]'^)“^^^ be approximated by 
reciprocal effective radial wavenumbers This might be so when excess curvature at low degree n is 
driven by narrower scale fluid motions, perhaps an ensemble of ephemeral eddies underlying the laminar 
sub-layer, that depend but weakly on multipole power at degee n. 

Similarly, for any set of PDFs assigned to PDFs assigned to 

^t ^n^^» define the ratio of expectation values ^ ^n^* characteristic time-constants 

In = nRncJ/fFncJ)^^ = ^ (14b) 

remain independent of radius in a source-free exterior. For example, if PDFs assigned to (gn*^. hj^”') were 
independent of m, and, in marked contrast to (8), SV were attributed to uniform westward drift at azimuthal 
velocity -WjjSinG alone, then {F^^} would be (nWo/c)^{R„j,}/3 and would be 3c^(nwQ)'2; however, all 
the R„, Fjj, and g^ll would be steady. 

3.3 Constant Aspect Ratio Hypothesis 

Effective radial wavenumbers k^^ for the core-source field can be compared with known horizontal 
wavenumbers To do so, introduce the square aspect ratio = k^^^/k^n^ and rewrite (14b) as 

^n ^ (ll^nc»/(Fnc))^^^ = fioc2[An n(n+l)]-l . (15) 

With no kinematical reason to prefer aspect ratios that either increase or decrease with n, constant = A 
seems the natural hypothesis to be tested. 


(L 


Some insight into this hypothesis is obtained by comparing radial diffusivity which is 

molecular (tio)'^ in the anisotropic laminar sub-layer, with lateral empirical diffusiviiies " 

Aj^/|XO. (These are "empirical" only in that they might eventually be determined from the ratio of geologic 
time-averaged observational spectra t^). The are amplification factors converting molecular diffusivity 
into lateral empirical diffusivities. They are also indicies of anisotropy. If there were but one lateral 
empirical diffusivity equal to a single lateral eddy diffusivity UL, then the single squared aspect ratio A 
would equal the magnetic Reynolds number paUL. This describes a scale invariant diffusive anisotropy. 
In contrast, the uniform westward drift example above has but one lateral eddy, yet many lateral empirical 
diffusivities, and would have be |iawQc(3)‘^^^(n+l)"^. 

More generally, many superimposed eddies of many lateral scales khi'^» speeds Uj(t), and circulation 
times I I can reconfigure flux from main field degrees j to other degrees n via (9). From a purely 

kinematic viewpoint, the many lateral eddy diffusivities I Uj(t)/kj^^ I do not explictly depend upon this S V 
degreen. Even an implicit dependence via selection rules (e.g., n~ I i +j| ) becomes obscured over time 
as eddies mix and remix all magnetic modes via interactions with intermediating modes. Moreover, an 
ensemble mean lateral eddy diffusivity, averaged over a kinematically unbiased ensemble of ephemeral 
eddies, cannot depend on the S V the eddies happen to induce. When, on average, neither radial molecular 
nor lateral eddy diffusivities { I UiCO/k^il } depend on SV degree, they cannot be used to construct lateral 
empirical diffusivities that do. The hypothesis of one such quantity, hence constant aspect ratio, remains as 
the seemingly natural null-hypothesis. Of course, detailed kinematical assumption (e.g., pure westward 
drift) or dynamical argumentation (e.g., strong Lorentz feedback) as to which particular eddies and 
magnetic modes might eventually dominate can suggest other hypotheses. 

This hypothesis is necessarily restricted to horizontal scales exceeding the thickness of the sub-layer 
itself (kfjp"^ » dg). It is clearly intended for magnetic scales (n < N£, n « Nj-)) broader than those 
dominated by molecular diffusion, hence time scales long compared with sub-layer delay lime (e.g., » 

padn^). Moreover, so long as fluctuations in and Fj^ are about the same magnitude as the means, a 
single epoch value for a single spectral ratio (e.g., R6(ti)/p6(ti)) could easily differ from by a factor of 
two or more. Such natural variability ought not be mistaken for a multiplicity of aspect ratios. 

Observational Check. Reliable, undamped geomagnetic field models through degree and order 12 
or more are available, albeit not for the geologically long time intervals needed to determine 
{Fj^(t)}/{Rp(l)l. Observational values of Fj^(a,t)/Rj^(a,t) were, however, calculated for epochs 1960, 1970, 
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and 1980 from undamped Magsat model GSFC 9/80 (Langel et al., 1982) fifteen years ago in hopes of 
deducing a dispersion relation for SV. The constant aspect ratio hypothesis predicts 

= (A/^ocV[n(n+l)]2; (16) 

therefore, at each epoch we fitted InA’ + pin[n(n+l)] by least squares to the observational ln(Fjj/Rjj) for 
degrees 3 through 12. The three resulting estimates of P average to 1.957 + 0.156 (la). This agrees with 
the constant aspect ratio prediction p = 2. It does not agree with the uniform westward drift example above. 

The three estimates of InA’ indicate (A’)'^^^ is 2,639 years to within a factor of 1.94 (2c); therefore, 
very recent field behavior suggests = 2,640/[n(n+l)] years. The extrapolation of this relation to the 
dipole, Tj = 1320 yr, agrees with the observational values [Ri(t)/F|(t)]^^ from GSFC 9/80, which average 
to 1275 yr. The extrapolated range for the quadrupole, i 2 = 220-880 yr, exceeds observational 
[R 2 ( 0 /F 2 (t)]^^, which average to 126 yr, perhaps due to a rapid rebound from diminished R 2 . 

With this Tp relation and core radius c = Cp^ = 3.5 Mm (from section 4), the lateral empirical 
diffusivity is about 150 m^/s. With (14b), and supposing kj.^ t n/c to allow forced dipole decay, this Xp 
relation also offers an upper bound on core magnetic diffusivity. The bound, (pa)'^ < 30 m^/s, is uncertain 
by a factor of two, as is the lower bound on core conductivity, a > 2.7x10^ S/m, obtained by assuming p = 
Po* The bound is consistent with a high conductivity liquid metallic outer core. 

When emerging S V is traced to underlying fluid motion, the bound may be interpreted as a turbulent 
diffusivity. To reckon molecular magnetic diffusivity, divide 30 m^/s by the ejfective magnetic Reynolds 
number Rpj = . If a molecular conductivity is 5x 10^ S/m is assumed, then Rp^* is 19 and square 

aspect ratio A is about 92, each to within a factor of two. 

Generalized Stevenson-McLeod Relations. On the constant aspect ratio hypothesis, equations (16) 
and (8) at the top of the sub-layer are combined to yield 

~ ^n^l^nc^*")! “ [poc^/An(n+l)]^[Cn(n + l/2)(n + l)](cA)^"‘‘^ 


= K (n + 1/2) [n(n + 1)]-1 (c/r)2n+4 

(17a) 

•» K [n(n + (c/r)2"+4 

(17b) 

= K(n+ I/2)-J (c/r)2n+4 

(17c) 


for 1 < n < Ng and r > c. With K set to Kj^(a/c)^ and r set to a, (17c) reduces to McLeod’s rule (6); 
however, (17) is obtained from constant aspect ratio, instead of a form for Pp(co), and does not exclude 


dipole power. This allows more extensive testing. Failure of (17) at narrow spatial scales (n « Ng) may in 
part be due to failure of (16) as well as (8). 

Particular spectral forms (7) and (17a-c) are of general form 

|R„^.(r)) = Kq(n)(c/r)2n+4. (18) 

The various q(n) approach Stevenson’s (1983) relation at high n; they are approximately 2/(2n+l) and are 
thus but minor variations on McLeod’s rule when compared with extending it to the dipole and to 
geologically long time intervals^ The latter correspond to the hypothesis that, over such intervals, the core- 
source dipole field at the top of the core might be exceptional compared with other multipole fields in its 
nearly axial orientation, but not in its magnitude. 

3.4 Statistical Core Field Hypothesis 

With the expectation value of core multipole power {Rj,c) given by (17), the statistical core field 
hypothesis advanced here is that, over geologic time intervals, normalized core magnetic multipole power 
(2n+l)Rj^^/{Rj^^} at radii r > c is distributed as chi-squared (x^) with 2n+l degrees of freedom. These 
particular trial probability distribution functions (PDFs) for the Rj^^ are 

^[(2n+l)R^^{R„^n = P2n+l(r) = [2 r(n+l/2)] (x^) exp[-x2/2] , (19) 

where T is the gamma function, the are from (17), and 1 < n < Ng. By (19), standard tables giving 
the probability of obtaining x^^X^ also give the probability of finding Rnc^^^nc^ ~ x2/(2n+l). Note 
that (19) does not postulate PDFs for each individual Gauss coefficient; moreover, {R^^} from (17) is not a 
plain exponential and does not level off at any depth. 

Distribution (19) neither requires nor prohibits equipaititioning of multipole power among the various 
orders m within degree n. It holds if Gauss coefficients of degree n are random samples of a population 
with a zero mean Gaussian PDF of variance (Rp^(a)}/[(n+l)(2n+l)]. It can also hold when these 
coefficients are not normally distributed. Indeed, there are an infinite number of PDFs for individual 
coefficients of any particular degree that yield (19) (see Appendix A). Distribution (19) may thus 
summarize core processes that cause fluctuations about a geologic mean magnetic energy density of 
{R^^(c)}/2p per degree, but distribute such energy abnormally among the orders within that degree due to 
geometric effects. Earth rotation is an obvious source of such anisotropy. 


In a dynamo, kinetic energy maintains magnetic energy against Ohmic dissipation. The Coriolis 
force does no work; it does not change the kinetic energy of a fluid parcel. Rotational polarization of 
(geologically) turbulent core flow might thus produce magnetic anisotropy without much deviation from 
(17). Anisotropy due to motional inductive effects resembling the Coriolis vorticity effect should 
distinguish axial from equatorial dipoles. For example, consider surficially geostrophic flow in the 
uppermost main stream: ^•Vx(^^xv) « 0 for bulk angular velocity Q. Then fluid downwelling implies 
poleward flow. By (9), fluid downwelling draws in or attracts magnetic field line footpoints to form 
regions of strong radial field (core spots); the accompanying poleward flow implies poleward drifting core 
spots. This provides a flux partitioning mechanism for axial dipole formation, growth, and fluctuation 
(Voorhies 1991; 1992), 

This mechanism for inducing planetary magnetic anisotropy relies on the kinematic boundary 
condition indicated by a rigid mantle; it might also operate where a stably stratified fluid bounds a 
convecting conductor. Whether or not it indicates deviations from (17), such as persistent quadrupole 
diminution (Stevenson, 1983) or a shift of magnetic energy from quadrupolar to dipolar configurations, 
remains to be seen. Other mechanisms, perhaps involving a solid inner core or a thermo-chemically 
heterogeneous mantle, might also partition core surface magnetic energy unequally among the harmonic 
orders within a particular degree for geologic time intervals. Hypothesis (16), expectations (17), and 
distributions (19) are compatible with such anisotropy. 

Earth’s squared dipole moment m^m is (47ca^/jO.^^[R j(a,t)/2]. Granting it is mainly from the core, 
(19) predicts that 3m*m/{m*m}, which equals 3Rj^(a)/{R2^(a)}, will be distributed as with three 
degrees of freedom. The probability #? ( I X I )tl I X I of finding I x I in the interval [\% \ , I X I I X 1 1 is 
prtdx^/dlxlMlxl , so (19) predicts a Maxwellian distribution of absolute dipole moments I m I . The 
predicted the root mean square dipole moment is 

RMSDM = = W{Rl^.(a)/2}l%o ; (20a) 

the predicted mean absolute dipole moment is 

IdmI ^ {Iml} = 47ta3{[Ri^(a)/2]l^}/^o 

= 47ca^[4{Rj^,(a)l/37t]l/2/^^ = (8/37t) I'^RMSDM ; (20b) 

and the most probable absolute dipole moment is (2/3)^'^^RMSDM. The Maxwellian distribution for 
normalized absolute geomagnetic dipole is quite compatible with an anisotropic, dominantly axial, dipole. 
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In following sections, one or two parameters of spectral forms {R^q} are estimated by least squares 
fittings to elementary functions of observational By (19), normally distributed residuals to the Rj^ fitted 
ought not be expected; moreover, the magnitude of these residuals should be large compared with the 
uncertainty in R^^ indicated by geomagnetic field model covariance. This is because the natural variability 
of the core-source field allowed by (19) turns out to be much larger than formal uncertainty estimates in the 
low degree R^^ derived from Magsat data. The modest decrease in the natural variance of the R^cCc) from 
(19) as n increases from 1 to 12, together with the Gaussian factor in (19), suggest least squares estimates 
closely approximate maximum likelihood estimates. Nonetheless, each estimate does not yield rigorous 
covariance for parameters (K,c); stated uncertainties summarize variances caused by geomagnetic field 
model and truncation selection. 

4 MAGNETO-LOCATION OF EARTH’S CORE 

If (17) is correct, then the two parameter fit of {R^^^(a)} to the spectrum calculated from spherical 
harmonic models of Magsat data should give accurate magnetic estimates of the core radius, denoted Cj^, 
and amplitude K. Such estimates were computed by least squares fits of spectral form (18) to observational 
R^(a) data. Specifically, In[{Rj^^(a)}/q(n)], which equals [nln(Cj^/a)^ + InK + 41n(Cj^/a)], was fitted to 
ln[Rj^(a)/q(n)]. The R^(a) data fitted were calculated from the degree 13 model GSFC 12/83 (Langel & 
Estes 1985) and, to reduce aliasing of crustal fields, the degree 60 model M102189 (Cain et al. 1990). 

The top three rows of Table LI give the computed using the q(n) listed in the first column. The 
second and third columns list Cj^ computed from Rj^(a) for degrees 3 through 12; the fourth and fifth 
columns list computed from R|^(a) for degrees 1 through 12, The fourth row gives leveling radii c* 
obtained using plain exponential form (3); the c* computed from degrees 2 through 12 are shown in 
parentheses. Such omission of dipole power raises c* by about 102 km. 

The top row of Table I.l shows that core radii determined by McLeod’s rule agree very well with the 
seismologic core radius of 3480 km. The closest agreement comes from the fit of (17c) to the R^W values 
from M102189 for degree 3 through 12. This is illustrated in Figure I.l, which graphs Rjj(a) (diamonds 
connected by dashed line segments), the fit (solid curve with Cj^ = 3484.5 km and ~ 

5.3649X 10^^ nT^), and the extrapolation to degrees 2 and 1 (fine dashed curve). The vertical bars attached 
to the fit show the 80% likelihood range deduced from per 2n+l degrees of freedom (19). (These are 
not errorbars attached to the data, which would be very small). As expected, two of the ten points fitted lie 
just outside the 80% likelihood range (Rg and R^). The low degree extrapolation shows R | is greater, and 
R 2 less, than expected; yet both are within the 80% likelihood range. 


Table I.l shows McLeod’s rule (17c) as most accurate (top row) and variation (I7b) as almost as 
accurate (second row); even form (7) gives core radii at most 1.8% above Cg (third row). The first two 
values in the top row average to 3488.2 + 10.3 km (2a). The four values in the upper left quadrant average 
to 3489.9 + 8.6 km (2o). The last two columns show that inclusion of dipole and quadrupole powers raises 
Cj^ by at most 22.3 km. That the effect on c^^^ is so small (<0.64%) is perhaps no more surprising than the 
efficient relaxation of Rj and R 2 towards expectation values. The effect of the somewhat lighter weights 
PDFs (19) would assign to lower degree Rj^ can be approximated by averaging results from degrees 1-12 
with those from degree 3-12. 

The fourth row of Table I.l shows that c* from plain exponential form (3) systematically 
underestimate Cg by about 200 km. It is doubtful that discrepancies this large result from crustal 
contributions to the Rj^(a) fitted. Indeed, coestimation of parameters describing recently advanced power 
law forms for the crustal spectrum R^x(^) would slightly decrease, not increase, c* and Perhaps 
geophysically important deficiencies of (3) have been overlooked due to omission of R| from previous fits. 

These spectral core magneto-locations are more accurate and precise than those obtained using the 
frozen-flux approximation alone (Hide & Malin 1981; Voorhies & Benton 1982; Voorhies 1984). The 
latter methods rely heavily upon uncertain phase information in harmonic orders m and, more importantly, 
upon uncertain secular change information from either SV models or main field models at different epochs. 
The spectral method allows for magnetic diffusion as well as motional induction and relies only upon the 
Rp spectrum at a single epoch. These advantages reduce scatter. Some 44 frozen-flux core locations 
obtained from diverse field models at various truncations of expansion (1) average to 3506.2 + 300.9 km 
(la, Voorhies 1984, equation (3.21)). The twelve estimates of c^^ in Table I.l average to 3504.8 + 37.0 km 
(2a). This is but 0.71% more than the seismic value. 

The success of Hide’s method nontheless suggests that secular change induced by main-stream fluid 
motion emerges largely unaltered by a visco-magnetic boundary layer. The success of the spectral method 
further suggests pseudo-random lateral magnetic transport by the main stream below, and diffusion of S V 
signal out the top of, this thin laminar sub-layer. Taken together, these inductive geomagnetic data analyses 
indicate a low viscosity, electrically conducting fluid core of 3.5 Mm radius without requiring seismic, tide 
and nutation, gravity, or astro-geochemical data to infer core composition, liquidity, and radius. 
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APPENDIX A: CHI-SQUARED FROM ABNORMAL DISTRIBUTIONS 
If 2n+l independent variables are drawn at random from identical, zero mean Gaussian probability 
distribution functions of unit variance, then the probability distribution function (PDF) for the sum of their 
squares, is well-known to be chi-squared with 2n+l degrees of freedom. It is evidently less well 

known that the reverse is not always true. For example, consider three independent real variables (X, Y, Z) 
on the open interval (-<»,«>) with PDFs 

x2exp(-x2/2) (Ala) 

^d(Y) = 5(Y) • (Alb) 

#?d(Z) = 5(Z) (Ale) 

where P is the bi-Maxwellian and 8 is the Dirac delta function. There is no chance of Y and Z being 
anything but zero, X can be either positive or negative, and + Y^ + Z^ = is distributed as 

= P(x2)d(x2) = 2 p^(X) l-^l d(x2) 

d(x2) 

= (27t)-l^ (x2)l/2 exp(-x2/2) d(x^) 

= [23/2 r(3/2)]-l (x^) 1/2 exp(-x2/2)d(x2) (Aid) 

where F is the gamma function. PDF (Aid) is chi-squared with three degrees of freedom. 

The difference between example (A I) and the usual case of three normal distributions with equal 
variance is important. With suitably normalized dipole coefficients (X g j^/D, Y ^ 0, Z ^ 0), (Al) 
describes a dipole with a zero mean bi-modally distributed axial component, variance {(gj^)^} = 3D^, and 
negligible tilt. This approximates a planetary magnetic field dominated by an axial dipole. The usual case 
(with (X, Y, Z) -> (gj^/o, g| Vo, hj Vo)) describes a dipole with no preferred direction at all and seems 
even less relevant to Earth than (Al). The isotropic case, with typical dipole tilts of order tan' V[2] ^^) = 
54.7"^, may be of interest for Uranus and Neptune; case (Al) may be of interest for Saturn (Connemey, Ness 
& Acuna, 1982); intermediate distributions are potentially more Earth-like and may be of interest for 
Jupiter (Connemey & Acuna, 1982). 
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There are an infinite number of sets of three PDFs for three independent variables which give the 
distribution with three degrees of freedom. To see this, consider the possibly singular distributions for 


independent variables (x, y, z) on the interval (-<», +») 


^ \ A I I 

pj(x) = A^lxl e 

I I '^y -y^/2 

#?2(y) = Aylyl ^e ^ 

, . A I I -z^/2 

#>3<z) = A^lzl e 


(A2a) 

(A2b) 

(A2c) 


where (A^^, Ay, Aj,) are normalization constants and (Pj^, Py, p^) are power law indicies which are less than 
one and sum to zero. Clearly 


p j ’(x^)d(x^) = 2p(x)\ — 1 d(x2) 

d(x2) 


= A^’lx 




(A2d) 


while similar PDFs for y^ or z^ follow by replacing x in (A2d) with y or z, respectively. The PDF for = 


2 _ 


x^ + y^ + z^ is obtained in the usual way via 


W 

= J i / #?l’(x^) ■ (x^+y^+z^)]<l(x^)d(y^)d(z^) 

0 0 0 


(A3a) 


and the fact that the offset delta-function is the inverse Lapalace transform (denoted La"^) of the 
exponential of its offset 

^2 . (x2+y2^_2)2 ^2 _s(x2+y2+z2) 

5[^2.(x2+y2+z2)] = La’ne ] = — J e e c 


2n\ 


ds 


(A3b) 


-loo 


where s is the Laplace transform domain variable. Substituting (A2d) and (A3b) into (A3a), and 
introducing u = v = and w = z^, yields 


AAA s^2 -d-i-p )/2 -u/2 

p(^2) ^ e La[u e ] 

2jci 

-(l+Py)/2 -v/2 -(l+Pz)/2 -w/2 

La[v e ] La[w e ]ds . 

With k„ = (1 - p„)/2 and 


(A4a) 


-(l-Hp„)/2 -u/2 k -1 -u/2 -k^ 

La[u e ] = La(u e ) = r(k^)[s + 1/2] 


(A4b) 
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it follows that 


A^’A^’A/ .i=o s^2 r(k^) r(ky) r(k^) 


^(^2) . _JLJL^ j 


2"' (s + l/ 2 )‘''''^V‘'z 


ds . 


(A4c) 


Provided k^, ky, and k^ are all positive definite and k^^ -i- ky + k^. = 3/2, 
#7(^2) = A^’Ay’A^’ r(k^) r(ky) nk^) La-l[(s + l/2)-3/2] 


r(k^) r(k ) r(k^) m 

= A^’A ’A^’ i (^2) e 

^ r(3/2) 

-^2/2 

= [23/2r(3/2)]-^^2)l/2g 


(A4d) 


Comparison of (A4d) with (Aid) shows that ^2 is in fact distributed as %2 with three degrees of freedom. 
The conditions that (k^, ky, k^) are all positive definite and sum to 3/2 are equivalent to the condition that 
(Px» Py’ Pz^ ^ There are an infinite number of triplets for (p^^, Py, p^) that 

satisfy these conditions. The form of PDFs (A2a-c) is A I x I 'P exp(-x^)] and is, in retrospect, one half the 
chi-squared distribution with typically fractional degree of freedom (1 - p) reflected about a zero mean; for 
p = (1 - 2k), the corresponding A is l/[2^r(k)]. 

It has been demonstrated that assigning the distribution chi-squared with 2n+l degrees of freedom 
to some variable does not necessarily imply that ^ is the sum of 2n+l normally distributed variables of 
equal variance. It particular, there are an infinite number of sets of three PDFs for the three degree 1 Gauss 
coefficients that imply the distribution of chi-squared with three degrees of freedom for normalized dipole 
power 3R|^{R } . Explict contact between (A2a-c) and PDFs for suitably normalized Gauss coefficients 

is omitted for brevity. Variances of PDFs (A2a-c) are restricted by the selection of power law indicies and 
are 2^^^r(k + l/2)/T{k). PDFs (A2a-c) also have singularities for 1 > py = p^. > 0. Such PDFs may be of 
interest in wave-particle physics, but are not the intermediate cases sought for describing planetary 
magnetism. 

Instead consider normal distributions for gj ^ and h j Then (Ala-c) makes it intuitively obvious 
that the PDF for g|^ yielding chi-squared with three degrees of freedom for the distribution of normalized 
dipole power will be a linear combination of Gaussian and bi-Maxwellian distributions. This is in fact the 
case, as shown below. 
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Denote [(a/c)^gi®, (a/c)3gj ^ (a/c)^h| h by [x, y. z] and the variances of (x, y, z) by (o^^, Oy^, 
o^'). The axial dipole hypothesis would have the means of y and z be zero, while isotropy of the equatorial 
dipole would have Oy and be equal. So replace (A2a-c) with distributions 

o 2-a 2 

p_(x)dx = (3/27tv2)l/2(3/v2)[-iL_^ x2 + a 2] exp(-3x2/2v2) dx (A5a) 

^ y2 J 

Py(y)dy = (27tay2)-l/2 exp(-y2/20y2)dy (A5b) 

p 2 (z)dz = (27ta22)'^'^ exp(-z2/2o22)dz (A5c) 

where v2 = aj^2 + cjy2 + 022. With ^2 s x2 + y^ + z2 = the problem is to derive (A5a) for p 

given normal distributions (A5b-c) and 

p(^2)d(^2) ^ (23/2r(3/2))-l[3^2/v2]l/2exp[-3^2/2v2](3/v2)d(^2) (A5d) 

subject to the constraint that the dipole is usually mainly axial: » Oy2 = For statstically 

independent (x, y, z) 


oo oo oo 

= I If P,,’(.yh - (x2+y2+z2)]d(x2)d(y2)d(z2) (A6a) 

0 0 0 ^ 

OO CX3 OO 

= (2710^0 )'U ! J p ’(x2)(y2z2)-l/2 exp(-y2/2o 2-z2/a 2) 

^ 0 0 0 ^ 

8[^2 . (x2+y2+z2)]d(x2)(i(y2)d(z2j (A6b) 

7 9 9 

or, using (A3b) and resetting u = v = w = 

ioo oo oo oo 

= {47^0^0 1 if! p^\u)txp(-su){\v/y^^exp{-sv -v/2a 

^ -ioo 0 0 0 ^ 

exp(-su - w/o^)exp(s^^) dudvdwds . (A6c) 

Recognizing and evaluating the Laplace transforms gives 

ioo 

P(^^) = (4Jt2ay02i)‘M La[pjj’(u)]La[v'^''2exp(-v/20y2)] 

-ioo 

La[w’^^exp{-w/2cy22)] exp(s^2)ds 

= (2naya2)-l La-l[La[Px’(u)] [r(l/2)]2[(s + l/20y2)(s + . (A6d) 


The Laplace transforms of (A6d) and (A5d) imply 


La[#;^’(u)] = ay02(3/2v2)l/2 [. 


(A7a) 


(s+l/20y2)(s+ 1/20,2) 

(s + 3/2 v2)3 


-]^/2(3/v2) 


In the special case Oy = 

P^’(x2) = o 2(3/2 v 2)I/2 La-1 ](3/v2) 

^ (s + 

and the inverse transform gives 


(A7b) 


o v2-3o 2 

(3A^2) [(x2)-l/2 + ^L(x2)-l/2]exp(-3x2/2v2) . (A7c) 

2r(3/2) Oy2y2 

Because x can be of either sign, p j^(x)dx is but half p ^’(x^)2xdx and (A7c) reduces to the linear 
combination of a zero mean bi-Maxwellian and a zero mean Gaussian (A5a). 

Although (A5a-c) are but one class of PDFs consistent with the distribution chi-squared with three 
degrees of freedom advanced for normalized dipole power in the text, they may be of some use in planetary 
and perhaps stellar magnetism. In (A5a) note that > 0; with the 

y A y Ay 

distribution enjoys the two peaks on either side of the mean (at gj® = 0) that correspond to normal and 
reversed axial dipole polarity, respectively. A suitable index of anisotropy, or ’tilt control parameter’, for 
PDFs (A5a-c) is e s [(c^^ - which is 1 for a purely axial dipole, zero for randomly oriented 

dipoles, and -1/2 for a purely equatorial dipole. For a purely axial dipole with o^ = o^ = 0, (A5a) would 
vanish at x = 0; then gj® could not pass through zero and reversals would be prohibited. So e « 1 describes 
inhibition of reversals. 

It is conjectured that PDFs (A5) describe the terrestrial dipole with fair accuracy. Very slow 
evolution of the variances Cy^ « c^) over geologic time is thought to reflect planetary evolution of 
the boundary conditions on the outer core (and vice-versa). PDFs (A5) may also describe other planetary 
dipoles (e.g., e is apparently about 1 for Saturn, somewhat less than 1 for Jupiter, and near zero for Uranus 
and Neptune). 
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Table 1. 1 : Estimates of Earth’s Core Radius in km from McLeod’s Rule and Other Forms 


q(n) 

GSFC 12/83 

M102189 

GSFC 12/83 

M102189 


3<n< 12 

3 < n < 12 

1 <n< 12 

l<n< 12 

[n + l/2]-l 

3491.8 

3484.5 

3501.2 

3496.0 


3493.3 

3486.0 

3506.8 

3501.6 

n'^ 

3511.4 

3504.1 

3543.1 

3537.8 

1 

3261.4 

3254.6 

3208.6 

3203.8 




(3311.8) 

(3306.0) 


Figure Caption 


Figure I.l shows the multipole powers Rj^(a) from model M102189 (Cain et al., 1990) for degrees 1 < n ^ 
12 (diamonds connected by dashed line segments), the two parameter fit of McLeod’s rule (16c) (solid 
curve with Cj^ = 3484.5 km and Kj^(a/cj^)^ = 5.3649x10^® nT^) to degrees 3 through 12, and the 
extrapolation to degrees 2 and 1 (fine dashed curve). Vertical bars show the 80% likelihood range from ^ 
with 2n+l degrees of freedom and are attached to the theory (solid curve), not the data. As expected, two 
of ten powers fitted lie outside this range. R j is greater, and R 2 less, than predicted, but both are within 


their 80% likelihood ranges. 






